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Abstract 

Using the relation of a set of nonlinear Langevin equations with 
reaction-diffusion processes, we note the existence of a maximal strength 
of the noise for the stochastic traveling wave solutions of these equa- 
tions. Its determination is obtained using the field-theoretical analy- 
sis of branching-annihilation random walks near the directed percola- 
tion transition. We study its consequence for the stochastic Fisher- 
Kolmogorov-Petrovsky-Piscounov equation. For the related Langevin 
equation modeling the Quantum Chromodynamic nonlinear evolution 
of the gluon density with rapidity, the physical maximal-noise limit 
may appear before the directed percolation transition, due to a shift 
in the traveling-wave speed. In this regime, an exact solution is known 
from a coalescence process. Universality and other open problems and 
applications are discussed in the outlook. 

1 Introduction 

The aim of the present note is to discuss some aspects of traveling wave solu- 
tions of a class of nonlinear equations with noise. These nonlinear Langevin 
equations are all related to the continuum limit of reaction-diffusion processes 
in statistical physics. 

*e-mail: robi.peschanski@cea.fr 
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A typical exemple of this class of equations is the stochastic Fisher- 
Kolmogorov-Petrovsky-Piscounov (sFKPP) equation. This equation reads 



dw(x, t) 



D 



+ A (w — w 



2 ) + e \Au(l - w) u(x, t) , 



(1) 



dt 



dx 2 



where the white noise verifies (u) = and 



(v(x, t)v(x', t')) = 6(x—x') 5(t — t') . 



(2) 



The deterministic part of ([I]) is the well-known FKPP equation pQ. D is 
the diffusion coefficient, A characterizes the activation strength since it is re- 
sponsible for an exponential grow of the FKPP solution in the dilute regime 
(when w 2 <C to). At some point the nonlinear damping term becomes ef- 
fective when entering the dense regime. The fluctuation contribution which 
sensibly modifies the solutions as we shall discuss in length, is characterized 
by the value of e, the noise strength. 

The same type of equation appears in various fields of statistical physics 
and, recently, in the domain of Quantum Chromodynamics (QCD), the in- 
teraction theory of quarks and gluons, where it happens to model [2] the 
evolution of the gluon momenta in the wave-function of a hadron or a nu- 
cleus when the energy increases. 

The main property of the FKPP equation [3] is to admit "universal" 
traveling wave solutions. Indeed, one finds in a quite general way traveling- 
wave solutions of the form w(x—vt) where the speed v is independent of the 
initial conditions (provided they are sharp enough in x). The reason is that 
a critical speed is selected after some evolution in time. Qualitatively, this 
property may be considered as a critical non-equilibrium mechanism since it 
corresponds to a solution flowing from an unstable (w = 0) to a stable (w = 1) 
fixed point. It results from the "strain" between the exponential increase 
raised up by the initial conditions in the linear regime and the damping due to 
the nonlinear terms. Subasymptotic terms of the speed and the structure of 
the wave front in the so-called "leading edge" region [lj can also be predicted. 
There are many results and applications of this traveling wave property in 
the literature [I]. 

In physical problems, one has often to take into account the existence of 
fluctuations. One simple exemple is the application of an evolution equation 
for a discrete set of "walkers" on a lattice generating statistical fluctuations 
at the edge, whose first approximation in the continuum limit is a cut-off, 
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e.g. on the function w of Eq.flTJ for the sFKPP equation. Indeed, it has been 
shown that, even with only a tiny cut-off [5], the effect on the asymptotic 
solutions is sizeable. The effect of fluctuations in the dilute part of the system 
(near w — 0) has been confirmed by simulations [6 J . It is amplified for two 
reasons. If of negative sign, a random fluctuation can drive the system to 
w = where it stays indefinitely, and thus acting as an effective cut-off. If 
positive, the fluctuation is amplified by the exponential increase in the dilute 
region and extends its effect far into the dense region. Both effects have 
been analyzed and clearly identified in the weak-noise regime [TJ. In fact, the 
solution may be interpreted as a stochastic combination of traveling waves 
with some dispersion A and the average solution changes from a function of 
x — vt to a function of In particle physics it corresponds to a transition 

from "geometrical scaling" [8] to "diffusive scaling" [9] properties of the gluon 
momentum distribution. 

In the strong-noise regime, the fluctuations are expected to have a decisive 
effect in low- dimensional problems. Indeed in the case of the sFKPP equa- 
tion, an exact solution can be found using a duality relation with reversible 
reaction processes at the coalescence limit [TUJITT]. It has also been applied to 
the QCD case [12] . In both cases, the exact solution features noisy traveling 
waves and correlations fully dominated by the strongest fluctuations. 

However, one may question whether one can go to arbitrary strength of 
the noise. Indeed, as we shall see through the general relation of the Langevin 
equations with react ion- diffusion processes, one has to take into account the 
existence of a phase transition which belongs to the universality class of the 
directed percolation (dp). At some point, the system no more propagates. 

Hence we are led for our purpose to take into account the fluctuations 
in an appropriate framework. The study of fluctuations is where the field- 
theoretical formalism is certainly useful, since they are properly taken into 
account in the very definition of the action. Our aim is thus to investigate 
what new information the field-theoretical formalism of reaction-diffusion 
processes can bring to the traveling-wave problem for the class of Langevin 
equations to which sFKPP and the QCD evolution equations belong. The 
field-theoretical formalism for reaction-diffusion processes [T3| [HI IS] has 
been already used in the context of the weak-noise regime of traveling waves 
|16j . where it confirmed the results of the cut-off analysis of Ref. [5]. Our 
aim in the present note is to use the field-theoretical formalism in the full 
range of noise strength, and more precisely at the strong- noise limit, where 
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renormalization effects allow for a proper account of the large fluctuation 
effects. 

In fact the main tool used of our note is to take benefit of a thorough field- 
theoretical analysis [H] of specific reaction diffusion process, the Branching 
Annihilation Random Walks (BARW), which we will relate to the strong 
noise problem of the class of Langevin equation we want to investigate. In- 
deed, the continuum limit of BARW processes leads to a Langevin equation 
of the same class. Hence the renormalization analysis of the BARW process 
leads to the existence of a maximal noise strength and to its evaluation. We 
will argue that this result extends to the class of langevin equations we study. 
However, the specific case of QCD will need a supplementary analysis, due 
to a different term in the corresponding Langevin equation. 

The plan of our note is the following: In section 2, we briefly recall some 
aspects of the Langevin equations applied to reaction-diffusion and QCD 
processes and in section 3, we infer the existence of the maximum noise 
strength from the study of the BARW process and give an evaluation. In 
section 4, we apply our result to the two exemples of the sFKPP equation 
and to the QCD nonlinear evolution. We close this note in section 5 by a 
summary of our results, with a discussion on the universality range of our 
result (for more general reaction-diffusion) and on possible applications to 
QCD. 

2 Langevin equations and reaction-diffusion 

processes 

Field theoretical formalism of the Langevin equations 

It is known that Langevin equations of sFKPP type can be analyzed in 
the form of a bosonic quantum field theory [13]. This formalism is particu- 
larly convenient to treat fluctuations superimposed to mean-field equations of 
FKPP type, since they are intrinsically taken into account by specific terms 
in the field-theoretical action. For instance, the field theory corresponding to 
Eq. (pQ) is defined in terms of bosonic fields w and their auxiliary counterparts 
w by the action 

S[w,w]= fd d x dt Iw [d t + DV 2 ] w-X ww(l-w)-— w 2 w(l-w)\ . (3) 
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In fact these auxiliary fields appear as external source fields for the deter- 
ministic (for linear terms in w) and the noise terms (for quadratic terms in 
w). Let us illustrate the field-theoretical formulation [13] with the example 
of Q: performing the path integral over the terms linear in the auxiliary 
field w gives rise exactly to the mean-field terms which, in the case of ([1]) is 
nothing else than the FKPP equation. Performing the gaussian path integral 
of the quadratic part in w 2 allows one to include the fluctuation contribution 
and generates correspondingly the noise term in the sFKPP version ([I]) of 
the Langevin formalism. 

This field-theoretical framework casts a useful bridge between the meth- 
ods of quantum field theory and the solution of reaction diffusion processes 
and with the treatment of equations of the sFKPP type (Ref. [H] will be of 
particular interest for us, while for a comprehensive review, see [15]). Indeed, 
in such reaction-diffusion processes it is possible to discuss the probability 
of evolution from an "inactive" state to an "active" one, under the effect of 
branching (or splitting), diffusion, annihilation and merging (or coagulation) 
of random walkers. In the continuum limit, it contains in full generality^ the 
field-theoretical description of the Langevin equations of the type ([T]). We 
put this connection on a precise footing in the following sections. 

One major interest of the field-theoretical treatment of the reaction- 
diffusion system is the treatment of fluctuations which play an important role 
in low-dimensional systems, such in our d = 1 case. It appears in particular 
through the renormalization of the field theory. In fact due to fluctuations, 
a phase transition happens which is in the universality class of directed per- 
colation [TilfTB]. 

Langevin Equation for QCD evolution processes 

In the QCD problem the region near the "inactive" state corresponds to 
the boundary of the dilute region for the gluon density (or "transparency" 
where the S- matrix S = 1) while the "active" state corresponds to the bound- 
ary of the dense region (or "strong absorption" where the S-matrix S = 0), 
see FigJTJ and Ref. [T7] for a review. 

The QCD evolution process which we refer to is thus the rapidity depen- 
dence of the tranverse-momentum gluon distribution in a nucleon or nucleus 

1 It is even more general since some reaction-diffusion processes cannot be described by 
a consistent Langevin equation (e.g. generating an imaginary noise term). 
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Figure 1: QCD traveling waves, from Ref.[21]. 

The traveling waves are mean field solutions of QCD evolution equa- 
tions, translationally invariant in time and joining the unstable fixed 
point (transparency or dilute medium of gluons), to the stable fixed point 
(strong absorption or dense gluon medium). 

target jTTj which evolves from the transparency to the strong absorption 
regime when the rapidity {i.e. the logarithm of the energy) increases. This 
equation, at least within an approximate framework can be given a field- 
theoretical interpretation in terms of non-equilibrium processes. 

Indeed, a coarse-graining approximation [18] leads to a simplified descrip- 
tion based on the stochastic equation 



dyT(L, Y) = X (-d L ) T(L, Y) - T 2 (L) + a/ kq. 2 s T(L, Y) u(L, Y) . (4) 



where L = logfc 2 , k is the gluon transverse momentum, Y the rapidity in 
units of the strong coupling constant a s and x the fraction of total energy 
brought by the gluon. We will consider x as the leading order in a s Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) kernel [19] 




(5) 
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Modifications by higher order corrections will not be introduced at the present 
stage of our derivation but they certainly occur and deserve further study. 
The noise u(L, Y) satisfies (u) = and 

(v{L, Y)u(L', Y')) = - 8(Y-Y') S(L-L') . (6) 

7T 



The deterministic part of (j4j) is the Balitsky-Kovchegov (BK) equation [20J . 
Its solutions have been found to converge to traveling waves after some ra- 
pidity evolution [2Tj . This fact is particularly clear if one uses the diffusive 
approximation as we shall see now. It leads to a deterministic part of the 
equation identical to FKPP with a modified noise term. 

Indeed, expanding the BFKL kernel to second order around some j c , 

X{l) =Xc + x' c (i-ic) + \x'c{l~lc) 2 = A + A a + A 2l 2 , (7) 

The value of 7 C is obtained through the equation 7 C x'(lc) = xilc) which 
determines the critical speed of the traveling waves. In fact, it can be shown 
[2T] that 7 C plays the role of an anomalous dimension characterizing the wave 
front formed by the "blocking" mechanism described in introduction. 

Introducing the space (x) and time (t) variables by the redefinition [21] 



T 

t = a s Y } x = L — Aia s Y, u(x,t) = — , (8) 

Aq 

Eq.fllj) gets mapped, within the diffusive approximation onto the Langevin 
equation: 

du(x,t) . d 2 u . , , . 9 , 2a s , — ; , , . s 

' =A 2 — + A (u(x, t) - u 2 (x, t)) + —^/kA u(x,t) u(x, t) . 9 
at ox z 7r 

It is thus easy to realize that the first term A corresponds to branching, 
A 2 to the diffusion constant. The space derivative term with coefficient Ai 
in (JTj), is thus translated into a time-dependent shift on the space variable 
by the redefinition (jBJ) of the kinematic variables. This shift plays, as we 
shall see, an important role in the application of reaction-diffusion processes 
to QCD and to the maximal noise problem in particular. 

Hence, in the coarse-graining approximation of QCD evolution process, 
the deterministic part of (El) can be interpreted as gluon momentum diffu- 
sion, branching and merging while the noise v is a simplified description of 



7 



the gluon-number fluctuation effects, which in the diagrammatic QCD frame- 
work is associated to the so-called "Pomeron loops" . The strength of these 
fluctuations is proportional to a s while the coarse-grained fudge factor k may 
be large. 
Writing 

D = A 2 , \ = A , e = —*J^A a , (10) 

7T 

Eq.Q is similar to the sFKPP equation (CQ), except for the quadratic term 
in the noise which is here absent. 

In the Langevin formulation, the quadratic term — w 2 in the noise of (flj) 
cancels the fluctuations near the stable fixed point of the mean-field equation 
at w = 1 , while they are not constrained near u = 1 in (Q . In the following 
discussion however, we shall consider Eq.([H]) and ([T]) together with the other 
ones differing only by the quadratic term in the noise as belonging to the 
same class of equations. The restrictive condition will be that no unphysical 
fluctuations be generated in the range spanned by the mean-field solution, 
i.e. the overall noise factor should stay positive in that range, such as for 
Eqs.flHD). 

Once this condition satisfied, the fluctuations are in general expected to 
be small in the "dense system" and without strong effect while they have an 
important role at small u or w which is the region of the dilute system. We 
will focus on this region in our study, and thus conjecture that our result 
is valid for the whole abovementionned class of Langevin equations with 
different quadratic term in the noise contribution. We will come back to this 
point in the discussion of our results. 

In the following we will use the diffusive approximation of the QCD equa- 
tion. However, one should keep in mind that a complete treatment of the 
fluctuations beyond the coarse-grained and fixed coupling approximations is 
still an open problem, since higher-order QCD contributions may change the 
features of the solution^. 

2 The effect of introducing the running of the QCD coupling constant in the QCD 
evolution equations with noise has been recently studied either through a numerical mod- 
elization [55] or on the BK equation with noise [23]. It seems that the fluctuation effects 
are transferred to higher rapidities than with fixed coupling. 
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3 Branching- Annihilating Random Walks 



As already mentionned the existence of a critical activation parameter has 
been found in BARW [Hj together with an estimate of its value within a 
2— d expansion at first order. An improved value could be found with higher 
order corrections or within the exact renormalization group approach [24]. 
We postpone this evaluation to further work. 

For sake of generality, let us start by considering a reaction-diffusion 
process in dimension d, including the transition processes A — > A + A, with 
rate A, A + A — ► (0, A) with rates (fi, //) respectively, and Ai — > A i±1 (where 
i ± i is a short-hand notation for nearest neighbours in dimension d) with 
diffusion rate D. One writes the following action [T3] for the probability 
distribution s of the species A 

S[s,s] = J d d x dt{s(d t -DV 2 )s-X n ss-Xs 2 s+(2fi+fi')ss 2 + (fi+fi')s 2 s 2 } 

(11) 

where takes into account the extra process A — > which is anyway 
imposed by renormalization of the A — > A + A. A general study of the 
reaction-diffusion action (fTTT) in terms of 4 — ei expansion [15] is possible but 
the critical dimension being d c = 4, (with the noticable exception of the 
BARW process) the study we perform in dimension d = 1 is not trivial. We 
will discuss further this point in the last section. 

Let us then consider more specifically the BARW transition processes, 
which has the interesting special feature [Hj to possess a critical dimension 
d c = 2. The BARW processes are described by the action ffTTj) restricted to 
y! = 0, i.e. no coagulation process A + A —>■ A. Using the same formalism 
[151 [16], an d considering for the present being the bare action (A^ = A), the 
theory corresponds to the Langevin equation 

^ = DV 2 s + As - 2/i s 2 + a/2As - 2/js 2 rj(x, t) , (12) 

This Langevin equation is well-defined i.e. the noise factor is real for s < A///, 
while the "active state" , which is the maximum value of s in the mean-field 
approximation, is twice as small as the value A/2// where the noise is again 
zero. The noise is non zero at the stable fixed point s = 1 of the mean- 
field equation but only very exceptional fluctuations could cause positivity 
problems under the square root. As already mentionned, we expect this 
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feature not to be relevant for our study. The Eg. fll2j) belongs to the same 
class of Langevin equations as flTJ and flH]) that we have considered, and thus 
the results obtained with the BARW process will be expected to be valid for 
the other cases. Through the redefinition 

e = yfljl (13) 



2/i ' 

we indeed recover the familiar form of these equations, namely the same as 
(JTJ), up to the quadratic term in the noise being now e yfw — ^ v(x,t). 

In Ref. [H], the fluctations are shown to renormalize the activation param- 
eter A through the one- loop effect A — > A + A —>■ 0, implying a counter-term 
in the action ffT2]) corresponding to the new transition A — > 0. As one may 
infer, the renormalization contributes with a negative sign, showing that the 
fluctuations tend to stabilize the classically unstable fixed point. 

Let us now show that the minimum obtained for the activation term in 
Ref. 1 14] corresponds to a maximal noise and evaluate its value. 

After field-theoretical renormalization [TJ] using the 2 — d expansion at 
first order, one finds the renormalized value \r of the activation parameter 

A -M ^gw^ )- (14) 

This induces, when A/j = and d < 2 a minimum value 

Amm — 7^ f 77—^ (15) 



D V2tt/ 

below which there is no transition from the inactive to the active state. In 
terms of the directed percolation phase transition which is met at this point, 
the system no more percolates. 

From the point of view of the Langevin formalism of (1121) in d = 1 and 
its solutions, it means that for 

A m in = 77 (-/r) ( 16 ) 



D V27T, 

there is an extinction of the noisy traveling waves. 

Considering; now the dimensionless noise e = e/( J DA) 1/4 , as defined in [IT], 
and the parametric relations (|13|) . one finds: 

£ "(53^ " c - = {dx-) =v ^' (17) 
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In fact, the key point of this maximal noise property is that when the noise 
is strong enough, the inactive phase (s = 0) from being originally unstable 
in the mean-field or with noise weaker than the maximal one becomes stable 
when the noise strength reaches the value (fT7j) . This prevents the instability 
which leads to the noisy traveling waves going towards the stable fixed point. 
This stabilisation by the noise is a phenomenon appearing under various 
forms in statistical physics (see e.g. [25]). 

Note that the value (jTTl) has been obtained from the evaluation [TT] of 
the minimum activation parameter at first order in the 2 — d expansion of 
the renormalization procedure. A more precise estimate could be probably 
obtained from higher-order calculations or from the exact renormalization 
group analysis |24J. The best hope would be an exact solution in al = 1. 



4.1 Maximal noise for sFKPP traveling waves 

In Ref. [IT], the sFKPP equation is studied by using duality properties with 
the A <-> A + A reaction-diffusion process with creation or recombinaison, 
each with its own rate. At strong coupling, this dual reaction-diffusion pro- 
cess goes to a limiting coalescence regime, where at most one walker occupies 
a lattice site and with a probability of giving birth to another one in a neigh- 
bouring site. This limit allows one to find exact analytic solutions [TU]. The 
resulting average wave speed, in terms of v c = 2a/ DA the speed in the ab- 
sence of noise (5j [TO] and e = e/(DX)i the dimensionless noise strength ffTTl) 
takes the form [TT] . 



The comparison with numerical simulations has been done in Ref. [TT] and 
reproduced in Fig.2. It shows the following features: in the range of weak 
noise, the normalized speed agrees quite well with the theoretical predictions 
[5]. With larger noise strength, one can observe that, when 0(1) < e < 
(9(5— 10), the numerically measured speed agrees reasonably well, within 
some statistical spread, with the theoretical prediction using the coalescence 
approximation. 

However the end of the spectrum, not discussed in Ref. [TT], seems to 
occur, compatible with a maximal noise limit. Indeed, the average speed of 



4 Applications to sFKPP and QCD 




(18) 
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the numerical solution appears to stop near, indeed slightly superior, to the 
2— d expansion value \J&k ~ 5, see Fig. 2. We thus note that the approximate 
value of e max is within the range of our estimate from the BARW process. It 
is an interesting challenge bpth from the numerical and theoretical point of 
view, to check whether more refined evaluations could confirm this finding. 
We hope to improve the this value by going beyond the 1-loop result or via 
the exact renormalization-group approach [23]. 

This region of noise strength below but near the transition where fluctua- 
tions stabilize the unstable mean-field fixed point corresponds to the critical 
region where the diffusion-reaction system meets the critical line of a phase 
transition of the directed percolation class [131 [23] . Hence, beyond the de- 
termination of the end-point of the noise strength, it would be interesting to 
discuss the other predictions coming from the renormalization properties at 
the phase transition [15] , for instance the relevance of the critical exponents 
to the description of the sFKPP solution and the matching with the near- 
by coalescence limit approximation of the related reaction-diffusion process. 
This study is certainly deserved for the future. 



4.2 Maximal noise for QCD traveling waves 

Coming back to the problem of QCD traveling waves with strong fluctuations, 
we have already noticed (see (jH])) that in the diffusive approximation of the 
QCD equation there is a change of space variable 

x = L-A x a s Y , (19) 

where A\ is the (negative) coefficient of the kernel expansion ([7]). 

Hence, in the absence of noise, Eq. fll9j) amounts to a shift for the wave 
speed 

nc •noise L ^^=] , (20) 



with respect to the critical value v c = 2^A A 2 . 

In the noisy regime dual to coalescence one finds [12 



stronq noise -, \ a \ , i 

V QCD .1 |^1 1 ^ 1 |A| 



< : 111 (21) 

v c (e) 2 2v^4^ " (imax) 2 2v/A^V V ' 

where we have given the limit corresponding to the maximal sFKPP noise. 
Following (I2T1) and the expression 017|) of the maximal dimensionless noise 
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strength, the normalized speed of the average QCD wave depends critically 
on the value of the ratio p = 2 ^a^a 2 com P are d to the maximal noise strength 

• If < p < 1 : The QCD speed goes to zero for a value of the noise 
strength which has not reached the value corresponding to ( fl7|) : The 
system does not reach the phase transition for positive speed (required 
by the physical interpretation of a rapidity evolving gluon distribution. 

• P < Knax '■ The system reaches the critical point where appears the 
directed percolation phase transition where the wave traveling stops. 
The behaviour of the solution near the endpoint should be given by the 
renormalization group analysis. 

In fact, if one sticks to the QCD leading-log order BFKL kernel ([5]) and the 
value Imax ~ 5 — 10 extracted from the simulations [IT] , one sees that we are 
in the case of the first item: the directed percolation transition point cannot 
be reached. Indeed, in Ref . [T2] . two options were discussed for the diffusive 
approximation of the kernel for the strong-noise region, depending on which 
value of 7 C the kernel has to be developped up to the second derivative. If 
one considers the same critical value 7 C as in the no noise case, one finds 
with (A Q ,Ai,A 2 ) = (9.55,-25.6,24.3) and thus p ~ .9 . If one considers 
a small value of 7 C 1 , one finds (A , A%, A 2 ) ~ (37^*, — 37~ 2 , 17c -3 ) and 
thus p = y/3/2, independently of j c . In both cases we find p « .8 ~ .9 < 
tmax ~ 5—10, and thus, referring to FigfJl the system of QCD traveling waves 
ends (reaching vqcd = 0) within the coalescence regime, without reaching 
the sFKPP maximal nois^E 

In order to give some appropriate estimate, using relation (flUI) . and writ- 
ing e < e^f, one finds the relation 

2 *V5<('^V. (22) 



7T V 

Hence, for ~ 0(1) the maximal noise could be reached even in the per- 
turbative regime (— small) if the fudge factor k is large enough. On the 
other hand, it is interesting to consider the effect of increasing the coupling 



3 The result could change with corrections from higher QCD pertubative orders but 
would need a profound modification of the kernel. One could also consider yet unknown 
strong a s coupling effects. 
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constant oc s . We see that in our framework it leads to an increase of the noise 
and consequently a decrease of the wave speed. Such an effect is otherwise ex- 
pected in QCD when considering the transition towards the non-perturbatice 
regime. Both considerations give an incentive to developping further a coa- 
lescence model starting from the original QCD formukation e.g. directly in 
terms of branching annihilation and merging of gluons or QCD dipoles. 

5 Summary of results and discussion 

In the present note, which used some known results on the field-theoretical 
treatment of stochastic Langevin equations of sFKPP type and of related 
react ion- diffusion systems, we have found the existence of a limiting maxi- 
mum strength of fluctuations, which can be quantified by a maximum value 
of the noise parameter in the Langevin equations. Below this value, there 
exists noisy traveling wave solutions which describe the transition of the sys- 
tem from an "inactive" phase to an "active" phase. Beyond this point, the 
fluctuations are strong enough to prevent the system to undergo such a tran- 
sition and it stays in the inactive phase, which is otherwise unstable. Hence, 
we have yet another striking example of random fluctuations contributing to 
stabilize a deterministic unstable point. 

Pursuing the analogies with the associated reaction-diffusion systems, the 
critical value of the noise corresponds to a directed percolation critical point, 
where the systems ceased to percolate. In this framework the noisy traveling 
waves could be interpreted as an average over the percolation path histories. 

We have discussed two applications of the results. In the sFKPP case, we 
interpreted the existing numerical simulations showing the compatibility with 
the existence of a maximal noise strength in the range predicted by the field- 
theoretical calculations. For traveling waves in the diffusive approximation 
of the QCD evolution equations, the traveling waves go to zero before the 
transition point of the sFKPP analoguous equation due to a shift of the speed 
(at least for the leading order kernel in a s ). The maximal noise-strength 
occurs in the region related by duality with a coalescence, which is accessible 
to exact noisy traveling wave solutions, as discussed in [T2]. For large enough 
QCD fudge factor, it may be reached even in the perturbative QCD regime. 
Otherwise, it may indicate the slowing down of rapidity evolution expected 
from the nonperturbative QCD regime. 

An important aspect of the discussion is the extent to which our result on 
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the maximal noise is "universal" , that is valid for a whole class of processes. 

In our study we have met different Langevin equations, namely Eq.flTJ 
corresponding to sFKPP, Eq. (Q corresponding to the coarse-graining of the 
QCD evolution equation in the diffusive approximation and finally Eq. (|12p 
corresponding to the BARW process and for which we could extract the value 
of the maximal noise. Since these Langevin equations are all well defined and 
differ only by a quadratic term in the noise which is not relevant for the region 
near the unstable mean-field fixed point, we have conjectured that the result 
on the maximal noise should be the same, and thus given by (jTTI) in the 
first-order approximation of the 2 — d expansion for the BARW process. 

Since each of these Langevin equations corresponds to the continuum limit 
of a specific reaction-diffusion process with different parameters for merging, 
splitting and annihilation, it is interesting to address directly the universality 
property in terms of a renormalization problem of the field theories for these 
three cases. By extension, one may investigate the validity of the obtained 
maximal noise strength for all reaction-diffusion processes pertaining to the 
universality class of the directed percolation processes [T5], for the class of 
processes for which the formulation in terms of a Langevin equation differs 
only by a quadratic term in the noise (unless leading to imaginary noise 
term), which is irrelevant near the directed percolation transition point. 

This task is not easy in the framework of the d c — d expansion since the 
critical dimension is d c = 4 in the directed percolation class, which is the 
universality class of the phase transition of these reaction-diffusion processes. 
It is a long way towards d— 1. However, a qualitative general argument can be 
given which relies on the topological property of random walks in dimension 
d < 2. Indeed, in that condition, any two random walks will have probability 
1 to meet again. Hence, if the combination of reaction-diffusion constants 
contribute to a large enough dimensionless noise strength, one expects a 
stabilization of the otherwise unstable fixed-point. Starting with the general 
expression (JTTJ), and following a similar discussion as for the BARW case, 
one is led to conjecture e max = \^8n + (higher order contributions). It is 
thus deserving to check this result in d=l@. 

4 Note added in proof; The special case d = has been reconsidered very recently. 
In Ref. [26] . the authors examine the diffusion-reaction systems with various couplings in 
the zero transverse dimension approximation. They find that only the reversible process 
A —> A + A, A + A — > A leads to an increase of the solution with time (or the amplitude 
with rapidity in QCD language). They conclude that only the sFKPP equation (in its 
reduction to zero transverse dimension) to which this equation is related gives a satisfactory 
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As an outlook one could propose two distinct directions for future work. 
In the sFKPP case, it would be useful to develop the connection between 
the regime just below the maximal noise with the results obtained by the 
renormalization group method near the directed percolation critical point. 
In particular the knowledge of the critical exponents should give interesting 
predictions for the behaviour of the averaged traveling wave and its eventual 
deformation at the critical point. From the QCD point of view, the exis- 
tence of a maximal noise-strength at least 10 time smaller than the sFKPP 
one, leads to interesting predictions from the coalescence process via duality 
properties [TUJ [EE]. It would be interesting to develop on the already known 
studies [12J by looking to dual properties similar to the coalescence limit 
more directly from the fromulation of QCD evolution processes in terms of 
merging, splitting and annihilation of gluons or QCD dipoles. 
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noise strength 



Figure 2: Wave speed as a function of the noise strength for the solutions of 
the sFKPP equation, from ref. YTT^ . 

Vertical axis: traveling wave speed normalized to the critical one 2\J DX of the 
deterministic FKPP equation. Horizontal axis: dimensionless noise strength 
e = e/(DX)i. Line on the left: the weak noise prediction [51 [16]; Line on 
the right: the strong noise prediction [10]. One observes a maximal speed 
around a noise strength of order 10. 
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